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Abstract 

An exactly separable version of the Bohr Hamiltonian is developed using a potential of 
the form u{(3) + m(7)//3 2 , with the Davidson potential u((3) = (3 2 + (3q/(3 2 (where (3 is the 
position of the minimum) and a stiff harmonic oscillator for 1^(7) centered at 7 = 0°. In 
the resulting solution, called exactly separable Davidson (ES-D), the ground state band, 
7 band and 0^ band are all treated on an equal footing. The bandheads, energy spacings 
within bands, and a number of interband and intraband B(E2) transition rates are well 
reproduced for almost all well-deformed rare earth and actinide nuclei using two parameters 
(/3o, 7 stiffness). Insights regarding the recently found correlation between 7 stiffness and 
the 7-bandhead energy, as well as the long standing problem of producing a level scheme 
with Interacting Boson Approximation SU(3) degeneracies from the Bohr Hamiltonian, are 
also obtained. 
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1 Introduction 

The Bohr Hamiltonian [lj has been at the foundation of the collective model description of 
nuclei for over fifty years. Numerous solutions have been proposed since its derivation by 
choosing different forms of the potential V(/3,7) and solving the corresponding eigenvalue 
equation either analytically or approximately. Recently, this approach has undergone re- 
newed interest, due in part, to the development of the concept of critical point symmetries 
(CPS). These models, E(5) [2] and X(5) [3], are special solutions of the Bohr Hamilto- 
nian designed to describe nuclei at the critical point of the shape/phase transition between 
vibrational and 7-soft or axially symmetric deformed structures, respectively. 

In E(5) [2], a 7-independent potential of the form u{(5) is used, leading to exact separa- 
tion of (3 from 7 and the Euler angles [I], while in X(5) [3J a potential of the form u{(3) +u( , y) 
is assumed, leading to an approximate separation of variables in the special case of 7 ~ 0°, 
achieved by ^(7) being a stiff harmonic oscillator centered at 7 = 0°. In both E(5) and X(5) 
an infinite square well potential is used as u(/3), in accordance with growing evidence from 
microscopic calculationsJSJEJE] that the potential at the transition point between different 
shapes should be flat. Model predictions for energy spectra and B(E2) transition rates are 
parameter free (up to overall scale factors) in E(5) [2J, while in X(5) the predictions related 
to the ground state band and the excited + bands are parameter free, but 7 bands contain 
the stiffness parameter of the 7 oscillator [3J, E] . 

It often happens that a successful, but simple, model or approach spawns new gener- 
ations of related approaches. This is especially the case if, despite its success, the data 
reveal certain, albeit perhaps small, discrepancies with the simple approach. A classic case 
of this is the simple formula for rotational spectra [T] , which led to a myriad of alternate 
formulas (see for example Refs. [9j [TOj [HI [12]), usually more and more parameters and, 
not surprisingly, working better. Of course, each such case ultimately entails a judgement 
as to whether the additional complications are worth the improved descriptions they yield. 

The case of critical point symmetries is no exception. Despite their simplicity (square 
wells in (3 along with flat or harmonic oscillator potentials in the 7 degree of freedom) and 
their success in describing transitional nuclei, it was immediately recognized that there were 
important discrepancies with the data as well. One, for example, occurs in X(5) where the 
predicted energy spacings in the excited O^ band are far too large [T3J [H] . 

Since the advent of these CPS, a number of alternate geometrical models have been 
proposed and their predictions worked out. Some of these share with X(5) an extreme 
economy of parameters, others have one additional parameter. These models can all be 
solved exactly, either analytically or numerically. Some are, in fact, essentially identical to 
the CPS but are solved exactly, while others involve alternate, presumably more realistic 
potentials. One example, which we shall refer to occasionally, is the so-called Confined 
Beta Soft (CBS) model [T5], US] which takes as its starting point from X(5) but allows the 
inner wall to move out to the radius of the outer wall. As the inner wall moves, the spectra 
change smoothly from X(5) to a pure rotor. Other potentials [TTJ [TS], [T9], which we shall 
not consider, utilize triaxial shapes with non-zero values for the minimum of the potential 
in 7. 

It is the purpose of this paper to explore a few of the most promising geometrical models, 
to compare their predictions with each other and with the data. We deal only with the 



axially symmetric case at present, that is, nuclei whose potentials in 7 are of harmonic 
oscillator type with a minimum at 7 = 0°. 

These models can be grouped into three classes: one, called X ex (5) is simply an exact 
numerical solution [20J of the Hamiltonian of X(5), without the approximate separation of 
(3 and 7 variables used in Ref. [3]. This type of exact solution has now become tractable, 
using the novel techniques introduced in Refs. [2H [221 E3]- It is worth mentioning that in 
Xex(5), the 7 stiffness parameter is involved in all bands, while in X(5) the ground and (3 
bands are independent of the 7 stiffness parameter . 

The other two classes each take advantage of a kind of potential that is exactly separable 
from the start. Such potentials have the form [3] 



ues(P,j) = u([3) + ^ (1) 
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where ES stands for exactly separable. The first of these uses the same u{(3) and 14(7) as 
X(5) itself - that is, the square well in (3 and a 7 dependent potential given by a harmonic 
oscillator in 7. This is the so-called ES-X(5) solution [23]. 

The second group of the exactly separable class of potentials uses the Davidson potential 
[25] in (3, namely 
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<(3) =P 2 + J 2 (2) 

where (3q is the free parameter and gives the position of the minimum of the potential 
in f3. The use of the Davidson potential with an approximate separation of variables has 
been discussed in Refs. [261 ET] . In the present work, we examine the Davidson potential 
with an exactly separable potential, which we call exactly separable Davidson (ES-D). By 
including a harmonic oscillator potential in 7, analytic solutions can be derived in this 
form to describe well-deformed, axially symmetric nuclei. These Davidson potentials, along 
with the X(5) potential, are illustrated in Fig. 1(a) for both the approximate separation of 
variables (left) and the exactly separable cases (right). In Fig. 1(b), the Davidson potential 
in just the (3 degree of freedom is illustrated for a few values of the (3q parameter. 

Before proceeding to a detailed discussion of the present solution of the Bohr Hamilto- 
nian, it is useful to put the present work in a context of other solutions to the Bohr Hamilto- 
nian. This Hamiltonian has been solved analytically in the 7-unstable case [u(/3, 7) = u((3)] 
using the Davidson potential of Eq. (2) as the /3-potential [28J , showing, that with increasing 
values of the f3$ parameter, a transition from the spherical vibrator to a rigid non-spherical 
7-unstable structure occurs. The link provided by 0(5) between the 7-unstable geometrical 
model and the 0(6) limit of the Interacting Boson Approximation (IBA) model [29J has 
also been previously studied [30]. Later, it was shown |31j that the above mentioned 7- 
unstable Bohr Hamiltonian with the Davidson potential is characterized by the symmetry 
SU(l,l)xSO(5), with SO(5) due to rotational invariance in the five-dimensional collective 
space, and with SU(1,1) due to the Davidson potential. If the potential is allowed to also 
depend on 7, no algebraic solution has been found, but it has been shown that numerical 
calculations converge much more rapidly in an SO (5) basis with (3$ 7^ than in the usual 
spherical basis with j3 = [2H [221 E3]: the relevant SO (5) spherical harmonics having 



been calculated in Ref. [22]. The correspondence between this approach, called the alge- 
braic collective model |23j , and the different limiting symmetries of the Interacting Boson 
Approximation (IBA) model [29] has been studied in Refs. [321 [33] . This powerful method 
has been recently extended [34] to the SU(l,l)xSO(N) case. 

In view of the above, the present work is an analytic, special solution of the Bohr 
Hamiltonian with a Davidson /^-potential appropriate for axially symmetric prolate de- 
formed nuclei (since the 7-potential is taken to possess a steep minimum at 7 = 0), while 
the earlier solutions of Refs. [2BJ ED] refer to a Davidson /3-potential in a 7-unstable frame- 
work. As a result, the present solution will turn out to be appropriate for the description 
of well-deformed axially symmetric nuclei, which comprise the bulk of well-deformed nu- 
clei, while the solution of Refs. [28l |3T] is appropriate for those 7-unstable nuclei between 
spherical and moderately deformed cases. 

There are several advantages in the present ES-D solution which we will consider in 
detail. As mentioned above, no approximation is involved in the separation of variables. 
As a result, all bands (ground, 7, and 0) are treated on an equal footing depending on 
two parameters, the Davidson parameter (3 (which is the location of the minimum of the 
potential) and the stiffness c of the 7 oscillator. Finally, the j3 2 term in the potential solves 
the spacing problem in the (3 band that plagues the infinite square well solutions. Of course, 
with a minimum in 7 at 0° and a relatively steep potential in 7, the model is applicable 
only to axially deformed rotational nuclei. 

Despite this constraint, it will be shown that the present solution provides good results 
for the spectra and B(E2) transition rates of almost all well deformed rare earth and 
actinide nuclei. Furthermore, it provides insights regarding the recently found correlation 
[35] between 7 stiffness and the 7-bandhead energy, as well as the long standing problem of 
producing an Interacting Boson Approximation (IBA) SU(3) degenerate level scheme 
within the framework of the Bohr Hamiltonian. 



2 The ES-D model 

Our starting point is the original Bohr Hamiltonian [Tj 
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(3) 
where /3 and 7 are the usual collective coordinates, while Qk (k = 1, 2, 3) are the components 
of angular momentum in the intrinsic frame, and B is the mass parameter. 

We assume that the reduced potential, u = 2BV/h 2 , can be separated into two terms 
of the form 



«09 l7 )=„C9) + ^ (4) 



as in Refs. [U [TTJ> ESI EH] where the Schrodinger equation can then be separated exactly 
into two equations. 

For the potential in 7 we use a harmonic oscillator 



u( 7 ) = (3c) V (5) 

and u{0) is taken as the Davidson potential 
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«(0) = P 2 + % (6) 



where (3q denotes the position of the minimum of the potential. As described in Appendix 
I, the resulting energy eigenvalues are given by 



EnL = 2n + 1 + y L(L + ^ K " + - A + ffi + 3C(n, + 1), n = 0, 1, 2, . . . (7) 
For K = one has L = 0, 2, 4, ... , while for K ^ one obtains L = K, K + 1, K + 2, 

Bands occurring in this solution, characterized by (n, n 7 ), include the ground state band 
(0,0), the /?i-band (1,0), the 71-band (0, 1), and the first K = 4 band (0,2). The relative 
position of all levels depends on the parameters (3 and C. (C = 1c is used in order to keep 
equations similar to those in Refs. [3] and [21].) All bands are treated on equal footing 
[36], in analogy with the SU(3) limit of the Interacting Boson Model [29] . 

Details on the calculation of B(E2) transition strengths are described in Appendix II. 
We note that the u(j) potential used in the Bohr equation has to be periodic, because 
of coordinate symmetry constraints pp. In Ref. [20] both the proper periodic potential 
(1 — cos 37) and the approximate form j 2 , appropriate for small 7, have been used, yielding 
similar results. A more detailed study of this issue has been recently carried out [37], leading 
to the use of spheroidal or Mathieu functions. Periodic 7-potentials involving cos 37 have 
been used in an early solution involving a harmonic oscillator for the /3-potential [38] . 
as well as more recently in the framework of the algebraic collective model [21] 1321 [33] , 
where their treatment is tractable because cos 37 is, within a constant, an SO(5) spherical 
harmonic with v = 3 and L = (where v the seniority and L the angular momentum) [21J. 
The potential esc 2 37, which is the partner of the infinite well potential in supersymmetric 
quantum mechanics [39j, has also been used recently [40] in the Bohr Hamiltonian for 
triaxial nuclei. It is certainly of interest to examine the consequences of the use of periodic 
7 potentials in the present approach in subsequent work. 

In the present paper we are going to follow Ref. [8], normalizing AK = transitions 
to 2^ — > 0^, and AK = 2 transitions to 2+ — > Of. In this way normalization difficulties 
vanish. 

The spectrum and B(E2) transition strengths of ES-X(5) are described for completeness 
in Appendix III. 

3 Numerical results and comparison to experiment 

3.1 Energy ratios 

In Fig. 2(a), the R4/2 = E(4f)/E(2f) ratio as a function of the parameter C is shown for 
the ES-D solution (for a few values of the Davidson parameter (3 ) and ES-X(5), as well 



as for the exact numerical solution [Xex(5)] of Ref. [20]. The parameter C is connected 
to the parameter a of the exact numerical solution [20] through C = ^y/a. In Figs. 2(b) 
and (c), the ratios R / 2 = E(0^)/E(2f) and R 2 / 2 = E(2+)/E(2f), corresponding to the 
normalized f3 and 7 bandhead energies, respectively, are shown for the same solutions. The 
ES-D solution with /3 = corresponds to the ES-X(5)-/? 2 solution of Ref. [24J. 

From Fig. 2(a), it is clear that the ES-D and ES-X(5) solutions are appropriate mostly 
for well-deformed nuclei, while the exact numerical solution [20] is also applicable to less 
deformed nuclei (including the a = 200 case [C = 9.428] which gives results similar to the 
original X(5) model |3j). This difference is due to the f3 2 term in the potential u(/3) +u(j) / f3 2 
used in the exactly separable cases. Within the ES-D solution, the rotational limit of 
R4/2 = 10/3 is closely approached already for (3 = 4. 

As seen in Fig. 2(b), the normalized (3 bandhead energy, -R0/2, has a large dependence 
on the parameter (3q and shows less variation with the stiffness parameter C, particularly 
for large (3q values. This dependence is reversed for the normalized 7 bandhead energy, 
R2/2, which varies only slightly for different /3 values, but has a large dependence on the C 
parameter. As a result, the R2/2 an d -R0/2 lines cross in ES-D at values of C increasing with 
(3q. These lines also cross in the numerical solution [20], but they do not cross in ES-X(5). 
This point will be further discussed in the next subsection. 

Concerning the results of the exact numerical solution Xex(5) of Ref. [20] used for 
comparisons in this and in subsequent sections, our assignment of levels to a particular 
band follows the same as given in Ref. [20]. In particular, 2+ corresponds to 2 2 for 
a = — 450 and to 2^ for a = 500 — 1000, while 2^ corresponds to 2g for a = — 450 and 
to 2 2 for a = 500 — 1000. Similarly, 4+ corresponds to 4 2 for a = — 650 and to 4^ for 
a = 700 - 1000, while 4^ corresponds to 4^" for a = - 650 and to 4^ for a = 700 - 1000. 
These assignments are related to avoided crossings, as explained in Ref. [20] . 

3.2 Relative spacings within different bands and relative posi- 
tions of bandheads 

In Fig. 3(a) the energy ratio 

_ jg(2+) - Jg(0+) 

12/3 " — eW) {) 

is shown for the solutions under discussion. The ratio is exactly 1 in the case of ES-D, 
irrespective of the value of the Davidson parameter (3q. This is due to the oscillator term 
in the Davidson potential which gives equal rotational spacings in the ground state band 
and the (3 band. Thus, the same holds for the energy ratio 

_ E(4j) - E(2j) 
RaP ~ E(Af)-E(2ty (9) 

shown in Fig. 3(b). In Fig. 3(c) the energy ratio 

_ g(4+) - E{2+) 
^ " E(4f)-E(2i) (10) 
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is shown, while in Fig. 3(d) the energy ratio 

£(0+) - E{2+) 
Rp-t = R 0/2 - R2/2 = E(2 + ) ^ 

is given. Abrupt changes in the predictions of X ex (5) are due to the avoided crossings of 
(2+ 2+) and (4+, 4+). 

Experimental data for the energy ratios shown in Fig. 3 are exhibited in Fig. 4. Since 
the current solution is only applicable to well-deformed nuclei, the data included in Fig. 
4 is limited to A > 100 and R4/2 > 3.00. Comparing Figs. 3(a) and 4(a) we see that in 
terms of the energy ratio R213, which compares the level spacing within the (3 band to the 
level spacing within the ground state band, most nuclei exhibit a ratio slightly less than 
1.0. This feature is most closely reproduced by the ES-D solution which predicts a ratio 
of exactly 1.0, independent of the value of the Davidson parameter (3 Q . The predictions of 
the ES-X(5) solution for R 2 p are higher by 50% or more, and the predictions of X e x(5) are 
even higher. (This is the well known problem of overprediction of the spacing of the (3 band 
in the X(5) model by a factor close to two [T31 EH], which can be resolved by replacing the 
infinite well potential by a potential with linear sloped walls [41J.) The same is seen for the 
energy ratio R±p in Figs. 3(b) and 4(b). 

From Fig. 4(c), it is clear that the majority of the data for the energy ratio i? 7 , which 
compares the level spacing within the 7 band to the level spacing within the ground state 
band, is centered around values of 1.0. The predictions for R 1 from each of the solutions 
overlap and are consistent with the range observed in the data. The ES-D solution gives 
the largest range of predictions, since this is a more flexible model (2 parameters) compared 
with the single parameter X e x(5) and ES-X(5) solutions. Overall, all three solutions yield 
reasonable predictions for the 7-band spacings in deformed nuclei. 

The experimental energy ratio Rp 7 , which is related to the relative positioning of the 
j3 and 7 bandhead energies, exhibits a wide range of values spanning positive to negative, 
as shown in Fig. 4(d). As a result, we expect that the solutions exhibiting both positive 
and negative values for this ratio, namely the ES-D solution for not very high values of the 
Davidson parameter /3 and X e x(5), should better reproduce this feature. 

Summarizing the above observations, all three solutions under consideration, ES-D, ES- 
X(5) and X e x(5), are found to give reasonable predictions for the 7-band spacing, while 
the ES-D solution yields predictions which most closely reproduce the /3-band spacing of 
most deformed nuclei. The ES-D solution (for not very high values of the parameter (3q) 
and the X ex (5) solution appear to reproduce the relative positions of the /3 and 7 bandhead 
energies in a number of nuclei. Thus, the ES-D solution provides the flexibility to describe 
a wide range of observables (spacings within the (3 and 7 bands, relative position of the f3 
and 7 bandheads) with not very large values of the Davidson parameter j3q- 

3.3 B(E2) ratios 

Having examined the main features of the energy spectra, we turn now to the study of the 
characteristics of the B(E2) transition rates. As mentioned in Sec. 2, in order to avoid 
normalization problems, AK = transitions will be normalized to the 2f — > Of transition, 
while AK = 2 transitions will be normalized to the 2+ — > 0^ transition, as in Ref. [8]. 



We include in this comparison the predictions of the original X(5) solution as well as the 
U(5) and SU(3) limits of the IBA. X(5) predictions for ground — > ground, (3 — > (3, and 
(3 — > ground transitions are taken from Ref. j3], while X(5) predictions for 7 — > 7, 7 — > 
ground, and 7 — > /3 transitions are taken from Ref. [8]. SU(3) predictions for the ground— > 
ground and (3 — > (3 transitions are obtained with the standard quadrupole operator of the 
IBA [25], while predictions for (3 — > ground and 7 — > ground transitions are obtained with 
the extended quadrupole operator containing the extra term (d'xs + stx d)^ |29j . 

Intraband transitions within the ground state and /3 band are shown in Fig. 5. Within 
the ground state band, Fig. 5(a), the ES-D predictions lie in between X(5) and SU(3) 
for most values of C. Again, for /3 values of 4 and larger, the SU(3) limit is almost 
exactly achieved. For the transitions within the /3 band, shown in Figs. 5(b),(c), the ES-D 
predictions lie between U(5) and SU(3), approaching the latter with increasing values of 
/3b. In the case of the ratio B{E2; 2^ -* Oj) / B(E2; 2f -> 0|), the X(5) predictions lie 
below the SU(3) value, and not between the U(5) and SU(3) values, as might have been 
expected, since they are related to the transition between U(5) and SU(3). As with the 
energy spectra, the (3 band predictions again exhibit the largest differences between ES-D 
and X(5). 

Transitions from the 7 band are shown in Fig. 6. The ES-D predictions are consistently 
close to X(5) or intermediate between X(5) and the SU(3) limit. This is particularly true 
for the branching ratios from the 7 band given in Figs. 6(b), (c). 

The (3 band to ground band transitions are shown in Fig. 7. For the transition from 
the 0^ state to the ground state 2f, the predictions of ES-D are intermediate between X(5) 
and SU(3) for most values of C. The decay from the 2^ state to the ground state band, 
Figs. 7(b), (c), shows some variation between the models, but all are similar in magnitude. 
These small differences become more evident when branching ratios are considered, as in 
Fig. 7(d). For the ratio, B(E2; 2+ -> 4f) / B(E2; 2+ -» Of), the X(5) predictions are 
nearly an order of magnitude larger than the SU(3) ratio (also the Alaga ratio) of 2.6. The 
ES-D predictions are again intermediate between X(5) and SU(3). 

The 7 band to (3 band transitions are shown in Fig. 8. The predictions of ES-D for 
growing (3q approach X(5). 

In summary, in (almost) all cases the ES-D predictions lie in general between the X(5) 
and SU(3) predictions, with SU(3) already approached at (3q = 4. 

3.4 Fits to specific nuclei 

A search has been made to find nuclei for which the ground state, (3, and 7 bands (up to 
the point of backbending or upbending in each band) can be well reproduced by the ES-D 
solution. Since the ES-D solution is appropriate only for deformed nuclei, the search was 
constrained to nuclei with i? 4 / 2 > 3.00. Considering all such nuclei in the rare-earth and 
actinide regions, we find that almost all nuclei with a known 0^ and 2+ state can be well 
described in terms of energies by ES-D, as shown in Table 1. The quality measure 
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used for evaluating the rms fits performed, remains less than one in most cases. Out of 
the ~ 60 nuclei which meet the above criteria, there are only two cases where the ES-D 
solution does not provide a good description of all the bandheand energies, namely 152 Sm 
and 154 Gd. These exceptions are not surprising, since these nuclei are well described by the 
X(5) model, which uses a "flat-bottomed" potential in the j3 degree of freedom. The ES-D 
solution, on the other hand, incorporates a potential which is much "stiffer" in the f3 degree 
of freedom. Thus, discrepancies between the data and the ES-D solution are expected in 
transitional nuclei and, indeed, may be used to point to nuclei with flat potentials in the (3 
degree of freedom. 

Several B(E2) ratios obtained with ES-D using the same parameters as given in Table 

1 are shown in Table 2, which includes all nuclei of Table 1 for which nontrivial information 
on relevant B(E2)s is experimentally known [15]. More detailed level schemes for 156 Gd 
and 232 Th are shown in Fig. 9, as examples of the quality of the ES-D solution to reproduce 
detailed spectra. 

As seen in Table 2, the intraband B(E2) ratios within the ground state bands are 
reproduced quite well for a majority of the nuclei, despite the fact that B(E2) values 
have not been taken into account in the fitting procedure. Also, the theoretical 7 — > 
ground B(E2) ratios are in very good agreement with the experiment values. However, the 
theoretical 7 — > ground B(E2) strengths, when normalized to the 2f — > Of transition, are 
much lower than the experimental ones. This could be due to the normalization difficulties 
mentioned at the end of Sec. 2, which disappear if ratios of 7 — > ground transitions are 
used. Moreover, the theoretical interband (3 — > ground B(E2) values are consistently an 
order of magnitude higher than the experimental values. 

3.5 Bandheads 

The ability of the present model to reproduce the general experimental trends of Rq/ 2 = 
E(Q^)/E(2t) and of R 2/2 = E(2+)/E(2f) as a function of #4/2 = E(Af)/E(2f) is shown 
in Figs. 10(a), (b). Predictions of the ES-X(5), Confined Beta Soft (CBS) [13 [16] and 
Xex(5) solutions are also shown for comparison. From Eq. (7) it is clear that the energy 
levels of the ground state and (3 bands depend only on the parameter combination /3q + 3C, 
thus in Fig. 10(a) only one curve appears for ES-D, with j3q + 3C increasing from left to 
right. From the same equation it is also clear that the levels of the 7 band depend on the 
parameter combination /3q + 6C. As a result, different curves are obtained for ES-D in Fig. 
10(b) by fixing (3$ to different values and varying the C parameter. 

The predictions for Rq/ 2 a s a function of R4/2 are more less the same for the ES-D, 
ES-X(5) and CBS solutions. The X e x(5) predictions for R / 2 are slightly higher and above 
the overall trend of the data. 

As discussed previously, in the CBS solution and other X(5)-related solutions, the band- 
head energy of the 7 band depends on a free parameter. In the present exactly separable 
(ES) solutions, it is treated on an equal footing as the (3 bandhead energy. The plot of R 2 / 2 
vs. -R4/2 reveals that a large set of data corresponds to the ES-D region with fi between 

2 and 4. The same set is also described quite well by the ES-X(5) curve. Figure 10(b) 
also reveals that the predictions of the ES-D solution for the 7 bandhead energy are only in 
agreement with the data for R4/2 values larger than 3.0. This is again related to the present 



solution being applicable only to axially symmetric well-deformed nuclei, since when the 
parameter C becomes too small, the approximation of an axially symmetric potential is 
no longer valid. On the other hand, X ex (5) provides a better description of R2/2 for -R4/2 
values between 2.6 and 3.0. 

3.6 Gamma-stiffness 

In Ref. [35J a correlation has been found between the gamma stiffness of the potential and 
the ratio R2/2 — E(2+)/E(2f), with the gamma-stiffness increasing stronger than linearly 
as a function of i?2/2- In the present model, the gamma-stiffness coefficient (3c) 2 is shown 
as a function of R2/2 i n Fig- 11(a). It is evident that a stronger than linear increase is seen, 
which varies little with (3 Q , at least for reasonable values of the latter, as indicated from 
Table 1. The specific points corresponding to the rare earth and actinide nuclei of Table 1 
are shown in Fig. 11(b), exhibiting the same trend. 

A short discussion is now in place on the qualitative correspondence between the two 
parameters (fio, C) of the present solution and those of the usual two-parameter IBA-1 
Hamiltonian [i6l EE7] 



H((,x) = C 



:w)^-J-q*-q* 



(13) 



where hd = (ft • d, Q x = (s^d + efts) + x(d^dY 2 \ Nb is the number of valence bosons, and 
C is a scaling factor. The above Hamiltonian contains two parameters, ( and Xi with the 
parameter £ ranging from to 1, and the parameter \ ranging from to — \/7/2 = —1.32. 
The IBA dynamical symmetries are given by C, = 0, any x f° r U(5), £ = 1, x = ~V7/2 
for SU(3), and ( = 1, x = for 0(6). As remarked in Ref. [35], stiffness is proportional 
to the IBA parameter x- Thus, in the present case (3c) 2 roughly corresponds to \x\- On 
the other hand, we have already seen that increasing (3 leads to the SU(3) limit, thus /3 
is in qualitative correspondence to (. It should be emphasized, however, that while the 
IBA Hamiltonian of Eq. (13) can cover the whole region from U(5) (-R4/2 = 2) to SU(3) 
(-R4/2 = 3.33), the ES-D solution provides reasonable results only in the narrow region of 
_R 4 / 2 between 3.0 and 3.33. 

3.7 Occurrence of SU(3) degeneracy 

A long standing problem has been deriving from the Bohr Hamiltonian a spectrum similar 
to that of the SU(3) limit of the Interacting Boson Approximation (IBA) model [29]. The 
main features of the spectrum should be: 

a) The energy spacings among the 2 + , 4 + , 6 + , ...levels within the ground, f3 and 7 
bands should be identical. 

b) Furthermore, the 2 + , 4 + , 6 + , . . . levels of the j3 and 7 bands should be degenerate. 
In the present model, the spacings within the ground and (3 bands are identical, because 

of the oscillator term in the u(/3) potential, as already seen in subsection 3.2. It is therefore 
sufficient to examine the conditions under which the 2 + , 4 + , 6 + , . . .levels of the j3 and 7 
bands are degenerate. 
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From Eq. (7) it is trivial to see that the energy spacings in the /3 and 7 bands become 
equal for any L if C = 4/9 (since in this case the 2>C term in the (3 band is counterbalanced 
by the —K 2 /3 + 6C term in the 7 band, which has K = 2). However, this observation is 
of little physical significance, since the values of C appropriate for actual nuclei, appearing 
in Table 1, are considerably higher. 

Figs. 3(b), (c) indicate that in general the spacings within the 7 band are lower than 
the spacings within the f3 band by about 20% for most /3q and C values of interest. Thus 
within the present solution one can only hope to reproduce a situation with approximate 
degeneracy for the first few even levels of the f3 and 7 bands. 

From Eq. (7), the requirement E(2p) = E(2+) leads to the condition 9C 2 — 80C — 
16/3q — 356/9 = 0, while the requirement E(4t) = E(4^) leads to the condition 9C 2 — 
80C — 16/3q — 1028/9 = 0. Similar conditions occur from the requirement E(Lt) = E(L+) 
for higher L. These conditions can be approximately satisfied simultaneously only for very 
large values of /3q, which are outside the region of physical interest according to the /3q values 
appearing in Table 1, since too high (3 Q would result in too high a 0^" bandhead energy. 

In this way one is led to consider what happens for a fixed value of R0/2 = E(0t)/E(2i). 
In this case Eq. (7) easily leads to 

»-(%*- £)'-*-!■ (14) 

Thus for a given Rq/2 one can minimize with respect to (3q the rms deviation between the 
even levels of the j3 and 7 bands 



<J/3 i7 (L r 




\ L max /2 - 1 S I E ( 2 t) 7 J ' (15) 



the value of C obtained for each /3 from Eq. (15). Numerical results shown in Table 3 
indicate that a reasonable degree of degeneracy is obtained for L max = 10 and R / 2 > 15, 
which is of physical interest, since the Rq/2 values in Table 1 extend up to 27 . In Table 4, 
the results of the fit to 232 Th, corresponding to the values reported in Table 1 are given. 
In the case of 232 Th, which is very close to the Rq/2 = 15 case reported in Table 3, one 
can see that cr^ 7 (L max = 10) = 1.142, while a e ^(L mauX = 10) = 0.593. Therefore, although 
the overall fit is quite good, the degree of degeneracy obtained from theory is less than the 
one indicated by experiment. One could conclude that the present solution does contain 
parameter pairs which correspond to an approximate degeneracy of the low- lying even levels 
of the (3 and 7 bands, while at the same time the spacings within the j3 band are identical to 
the spacings within the ground band, however the problem of reproducing a SU(3) spectrum 
from the Bohr Hamiltonian remains conceptually open. 

3.8 Alhassid-Whelan arc of regularity 

It has been recently suggested that an experimental confirmation [IS] of the Alhassid- 
Whelan arc of regularity [49J, connecting the U(5) and SU(3) symmetries in the sym- 
metry triangle [50] of the Interacting Boson Approximation (IBA) model [29] is man- 
ifested in nuclei in which the j3 and 7 bandheads, 0^" and 2+, are nearly degenerate. 
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From Eq. (7) the requirement E(0t) = E(2+) leads to the condition 9C 2 — 68C — 
16/3q — 224/9 = 0. Given the fact that C has to be nonnegative, this condition leads 



to C — (34 + y 1380 + 144/3q)/9. Among the nuclei listed in Table 1, the ones satisfying 
the condition \E(2+) - E(0p\/E(2+) < 0.05 [48] are 158 Gd, 158 Dy, 170 Er, 178 Hf, 236 U, and 
248 Cm. 

From the /3q and C values listed in Table 1 , one can see that the above condition is closely 
fulfilled. However, since the ES-D solution is applicable mostly to nuclei with i? 4 / 2 > 3.0, 
the above mentioned condition describes only a small part of the arc of regularity close to 
the SU(3) limit. 

4 Conclusions 

In the present paper, an exactly separable version of the Bohr Hamiltonian, called ES-D, 
which uses a potential of the form u(/3) + u^/fi 2 , with a Davidson potential [3 2 + /3q/ ' (3 2 in 
the place of u(/3), and a harmonic oscillator with a minimum at 7 = 0° as ^(7), is developed. 
All bands (e.g., ground, j3 and 7) in this solution are treated on an equal footing, depending 
on two parameters, the Davidson parameter /?o and the stiffness c of the 7-potential. The 
solution is found to be applicable only to well deformed nuclei (with R±/ 2 > 3.0) due to the 
(3 2 denominator in the ^(7) term. Nevertheless, it reproduces very well the bandheads and 
energy spacings within bands of almost all rare earth and actinide nuclei, with R4/2 > 3.0, 
for which available data exists, as well as most of the inter-ground and intra-7 band B(E2) 
transition rates. The most glaring discrepancy concerns B(E2) values for the f3 band to 
ground band transitions which are typically overpredicted by an order of magnitude. The 
two exceptions where ES-D does not provide a good description of energy spectra are 152 Sm 
and 154 Gd, which have previously been shown to be well reproduced with the infinite square 
well potential of the critical point symmetry X(5). Furthermore, the ES-D solution provides 
insights regarding the recently found correlation between the 7 stiffness and the 7-bandhead 
energy, as well as the long standing problem of producing a level scheme with IBA SU(3) 
degeneracies within the framework of the Bohr Hamiltonian. 

However, several open questions remain, in particular, concerning the discrepancies in 
the B(E2) predictions. The underprediction of the 7 — > ground and 7 — > /3 B(E2)s can be 
attributed to two reasons. First, the /3 2 denominator in the 14(7) term, "pushes" the nucleus 
to more rigid axial behavior. This can be investigated through a detailed comparison of 
B(E2)s predicted by ES-X(5) and the exact numerical solution of Ref. [20j, since the same 
u(f3) and 14(7) potentials are used in both cases. Work in this direction is in progress. The 
second reason is the use of a harmonic oscillator potential for 14(7), as an approximation 
valid for small 7, instead of a potential periodic in 7. This can be studied through the use 
of a periodic 7 potential [37] in ES-D, since no approximations will be present in this case. 

Furthermore, an exact numerical solution parallel to Ref. [20] utilizing a u([3) +u( r y)//3 2 
potential with a Davidson potential as u{(3) should demonstrate the degree of importance 
of /?-7 coupling when compared to the present results. 
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Appendix I: Spectrum of ES-D 

One seeks [3] solutions of the relevant Schrodinger equation having the form \f(/3, 7, di) = 
4>k(Pi l)^M K^fidi where 9i (i = 1, 2, 3) are the Euler angles, V{6i) denote Wigner functions 
of them, L are the eigenvalues of angular momentum, while M and K are the eigenvalues 
of the projections of angular momentum on the laboratory-fixed z-axis and the body-fixed 
z'-axis respectively. 

As pointed out in Ref. [3], in the case in which the potential has a minimum around 
7 = one can write the angular momentum term of Eq. ([3]) in the form 



E 



Q 



fc=l,2,3 sin' 



( 7 -¥* 



M i + Q i + Qi) + Q -^- 

3 Vsin 7 



4 



(16) 



Using this result in the Schrodinger equation corresponding to the Hamiltonian of Eq. 
(jBJ), introducing [3] reduced energies e = 2BE/h 2 and reduced potentials u = 2BV/h 2 , 
and assuming that the reduced potential can be separated into two terms of the form 
w(/3, 7) = u(/3) + u(j)//3 2 , as in Refs. [U \TT\ [TBI HH], the Schrodinger equation can be 
separated into two equations 



1 d oi d L(L + 1) , , A 



UP) = e£ L (P), 



(17) 



1 



1 



d . d K 2 

^— -— sin 37— + — - I - T -^ r - 

sm 37 0*7 0*7 4 V sin 7 



4 N 



«(7) 



^(7) = ^(7). 



Eq. ( fT8l) for 7 « can be treated as in Ref. [3], considering a potential of the form 
u il) = (3c) 2 7 2 and expanding in powers of 7. Then Eq. (fl8|) takes the form 



1 9 d K 2 , o , 2 2 

ry I L (3c)^^/ Z 

7 $7 c?7 47 2 



^(7) = <W(7)> 



(19) 



with e 7 = A + ^-. The solution is given in terms of Laguerre polynomials [3] 



(3C)(n 7 + l), 



»i 



0, if = 0; 



ru 



K 



C = 2c, 

±2; 



»-, 



n. 



^ 7 ,|X|(7) = a 7 ,|x|7^ /2| e- (3ch2/2 4 X/2| (3c7 2 



= 0,1,2,..., (20) 

if = 0,±4; ..., (21) 

n = (n 7 - |if/2|)/2. (22) 



Eq. (TPTj) is then solved exactly for the case in which u(/3) is a Davidson potential 



«(/*) = ^ + f ■ 



(23) 



where /?o denotes the position of the minimum of the potential. In this case the eigenfunc- 

tions are [42] 

1/2 

Fm- -7 — ^ (3 a L a n + Hf3 2 )e-^ 2 , (24) 



V [n + a+ § 
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where T(n) stands for the T- function, L^(z) denotes the Laguerre polynomials, and 



2 V 3 1 

while the energy eigenvalues are 



3 '- /L(L+1) - A %? + « + 3CK + 1 ), (25) 



E n , L = 2n + a+^ = 2n + l + V l + | + Po + 3C(n 7 + 1), n = 0, 1, 2, . . . 

(26) 

For K = one has L = 0, 2, 4, ... , while for K ^ one obtains L = K, K + 1, K + 2, . . . 

In the above, n is the usual oscillator quantum number. A formal correspondence 

between the energy levels of X(5) and the present solution can be established through the 

relation 

n = s-l. (27) 

It should be remembered, however, that the origin of the two quantum numbers is different, 
s labelling the order of a zero of a Bessel function and n labelling the number of zeros of a 
Laguerre polynomial. In the present notation, the ground state band corresponds to n = 
(s — 1). For the energy states the notation E StL = E n+ljL of Ref. j3] will be kept. 
The full wave function reads 

*G0,7,0,) = FnW)rin„\K\(l)^ K m, (28) 

and should be properly symmetrized jl3] 



¥(0 l7 ,0«) = F^P) Vn ^ Kl (j)J - ]f + \ (V L M>K + (-1) L V L M> _ K ) . (29) 

V 167r(l + Ok,o) v ' 

It should be noticed at this point that Eq. (|T7|) for A = takes the form appearing in 
the framework of a X(5) solution with the infinite well potential replaced by a Davidson 
potential, called X(5)-D in the usual terminology. From the expression for A given below 
Eq. (TI9~1) it is clear that A = is achieved for K = and e 7 = 0, i.e. C = 0. It is therefore 
proved that the numerical results of the ES-D solution for the K = bands (ground state 
band and beta bands) will coincide with the corresponding results of X(5)-D. This result 
should be considered as a numerical coincidence, because C = is not acceptable in the 
framework of ES-D, since the approximation of 7 being close to zero collapses in this case. 
The lowest R4/2 value within the present model is obtained for /3q = and C — 0, which 
corresponds to the X(5)-/3 2 solution [51], giving i? 4 / 2 = 2.646. Thus, while /3o suggests a 
spherical shape, the contribution from centrifugal term in the potential results in a non-zero 
value for the average deformation. 

Appendix II: B(E2) values of ES-D 

B(E2) transition rates 

j^U.^.lgpg (30) 

v ' ; 16tt 2L + 1 v ; 
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can be calculated using the quadrupole operator [3] 

where t is a scale factor, and the Wigner-Eckart theorem in the form 



(31) 



(L'M'K'\T {E2) \LMK) = l (L2L'\MfiM')(L'K'\\T {E2) \\LK). (32) 

\/2L' + 1 

In ground — > ground, (3 — > ground, /?—>/? and 7 — > 7 transitions, only the first term of Eq. 
(j31j) contributes, since the relevant angular momentum coupling coefficients involving the 
second term vanish, while in 7 — > ground and 7 — > /3 transitions only the second term of 
Eq. (13TT) contributes, since the relevant angular momentum coupling coefficients involving 
the first term vanish. The final result reads 

B(E2; nLn.K -> n'L'n^K') = -^-t 2 ((L2L'\K, K> - K, K')) 2 B 2 n>L>nf>L ,C 2 nyiK ^ KI , (33) 

where 

B n ,L,n>,L> = J fiFf (/3)F% (f3)P 4 df3 (34) 

is the integral over (3, while C n K,n' ,k' is the integral over 7, in agreement to Ref. [8]. 
In ground — > ground, (3 — > ground, (3 — > (3 and 7—^7 transitions (AK = transitions), 
the integral over 7 becomes C ni) K,ri ,K' = fin-,,ri 3k,k', since (considering cos 7 ~ 1) it 
corresponds to the relevant orthonormality condition of the 7 wavef unctions, while in 7 — > 
ground and 7 — > (3 transitions (AK = 2 transitions) this integral has the form 

Cm,K,n^,K> = J Bin 77/^,1^1^,1^/1 1 sin 37^7, (35) 

since the volume element is pQ 

dr = (3 A I sin 37 1 sin 9dj3d'yd9d(j)dip- (36) 

For the bands considered here one needs the special cases of Eq. f[2"2"j) 

%,o = ^o,oe- (3ch2/2 , m,2 = C li27 e-( 3 ^ 2 / 2 , (37) 

where the Laguerre polynomials are unity since h = in both cases, as seen from Eq. (j2"2"l) . 
the relevant normalization conditions being 



(C ,o) 2 J e"( 3c ^ I sin3 7 |rf7 = 1, (Ci, 2 ) 2 J 7 2 e" (3ch | sm3 7 M7 = 1. (38) 

Then Eq. ( 135)) takes the form 

Ci.2,0,0 = Co,oC li2 |7 2 e -( 3c ^| sin3 7 |rf7, (39) 
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in which the integral is the same as the one appearing in the second normalization condition 
in Eq. (138|) . resulting in 

Ciao.0 = ^- (40) 

W,2 

Using the approximation | sin 37 1 ~ 3 |^y | and the integral 

r(*±i) 



x m e- ax2 dx 



2a 2 



the normalization conditions give 

(O ) ,o) 2 = 2c, ((7,,>)- 6r 



^ .-).. ,/-- ,2 ft. .2 Co.O 



C1.2 v^c 



The normalization is consistent with the one used by Bohr pQ. The same approximations 
are also used in Ref. 



Appendix III: ES-X(5) 

In the case of the ES-X(5) solution [24], in which u(/3) is an infinite well potential 

uW = l° "/>£Aj , (41) 

v ' \ 00 for (3 > (3 W 
the /3-equation becomes a Bessel equation with energy eigenvalues [3] 

ep, s ,L = {K L )\ k s , L = ^, (42) 

where x s> l is the s-th zero of the Bessel function J v {k s ^ L f3) with 



while the relevant eigenfunctions are 

ZsAP)=C s , L [3- 3/2 Uk s>L f3), (44) 

where C St L are normalization constants, determined from the condition 

rPw 

P%MdP = 1, (45) 





leading to 



a 2 r 2 ^+ 11 



y s,L 
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The full wave function reads 

*G0,7,0i) = C^p-^Mk^P^K^V^M), (47) 

and should be properly symmetrized 



¥(0,7, ft) = C^^Uk^)^^^)!^^ 



2L+l ^u K + (-l) L V"). (48) 



In calculating B(E2)s, the integrals over 7 and the Euler angles remain the same as in 
Appendix II, while the integrals over f3 take the form 

Bs^jj = C S , L C S ,, L , hj v {k S)L p)J v ,{k s ,, LI (3)(3d(3, (49) 



where the formal correspondence n = s — 1 holds. 

It should be noticed at this point that for C = the numerical results of ES-X(5) for 
the ground state band and the beta bands coincide with the results of X(5), as it can be 
seen from Eq. (133]) . As discussed at the end of Appendix I, this should be considered 
as a numerical coincidence, because C = is not allowed in the X(5) framework, since it 
destroys the 7 ~ approximation. 
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Figure 1: (a) Potentials in both the (3 and 7 degrees of freedom for X(5) (top) and the 
Davidson potential with (3$ = (middle) and f3$ = 2 (bottom). Potentials are shown for 
the approximate separation of variables (left) and the exact separation of the variables 
(right), (b) Davidson potential in the (3 degree of freedom for a few values of the parameter 

A). 
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Figure 2: (Color online) (a) The i? 4 / 2 = E(A~{) / E{2{) ratio as a function of the parameter 
C for the ES-D solution (for different values of the Davidson parameter /3 ), the ES-X(5) 
solution, and for the X ex (5) solution [20]. The X ex (5) parameter a is connected to param- 
eter C of the present solution through the relation C = (2/3) y^. (b) Same as (a), but 
for the ratio R0/2 = E(0t)/E(2f), corresponding to the normalized /3 bandhead energy, 
(c) Same as (a), but for the ratio _R 2 /2 = E(2^)/E(2i), corresponding to the normalized 7 
bandhead energy. See subsec. 3.1 for further discussion. 
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Figure 3: (Color online) Energy ratios (a) R 



2/i 



E{2+)-E(0+) 



fi(2+) 



£^[Eq. ©],(b)i^ 



[Eq. 



fc)iL 



E(4+)-E(2+) 



[Eq. ([TDD], and (d) R 



Pi 



E(0+)-E(2+) 



E(4+)-E(2+) 
£(4+)-£(2+) 



[Eq. (TTTj) ] as functions 



of the parameter C, for the same solutions shown in Fig. 2. See subsec. 3.2 for further 
discussion. 
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Figure 4: Experimental data for the same energy ratios shown in Fig. 3 . For each ratio, 
all nuclei with A > 100 and R4/2 > 3.0 for which sufficient experimental data (taken from 
Ref. [41)] ) exist, have been taken into account. The predictions for ES-D, ES-X(5) and 
Xex(5) are indicated in (a), (b), and (c). The X ex (5) predictions lie off scale to the right 
in (a) and (b). See subsec. 3.2 for further discussion. 



23 





2.0 


_ a ) 


i 








1 1 ' 1 


- 


1- — 

o 










-+— U(5) 


t 


1.9 












■ r n 


- 




P - u 


w 

PQ 


1.8 
1.7 












-•—(3 =2 

r 

-▼— P o = 4 

-°— X°(5) 
-*— SU(3) 


- 


T 

m 

PQ 


1.6 
1.5 




□— 

■ --. 

•- 


□— 

•— 


— n— 

— ■— 
— •- 


— D- 

— ■- 
— •- 


-n — n— n— n-n-n-n-n-n-n-n-n-n-n-n-n 


- 




1.4 




i 


M— 


— W— 


— W- 


till 






b) 


I 








' 1 ' 1 




+ — 
o 

T 

+ — 


1.5 


- 












- 


W 


1.2 


- 


■ 

•- 


-•— 




~m~ 
— •- 




- 


+ to. 

o 




■ 












- 


T 

+ ca 
(N 

CN 

w 
ffl 


0.9 


- 












- 


0.6 
2.0 


- 


1 


□— 


— n— 


— n- 


-n — □— □— n-n-n-n-n-n-n-n-n-n-n-n-n 

i 1 i 1 


- 


+"~ca- 
O 


c) 


1 








1 1 * 1 


- 


T 




■ 












■ 


<N 
CN 

w 


1.8 
















cq 


1.6 


_ 












_ 


+ ca 

<N 






□— 


□— 


— □- 


— □- 


— □ — □— □— □-□-□-□-□-□-□-□-□-□-□-□-□ 




T 

+ ca 


1.4 


- 


■ - 


=■= 

— i— 


=■= 
=w= 


=1= 


=1=1=1=1=1 iiiimiiii 


- 


pq 


1 i 


- 


1 








i . i 


- 



10 



15 



20 



Figure 5: (Color online) Intraband B(E2) ratios for the ground state and (3 bands vs. the 
parameter C as predicted by the ES-D model (labelled by the value of the (3q parameter), 
compared to U(5) [29], X(5) [31 E], and SU(3) [29] predictions, as described in subsec. 3.3. 



24 



+ -H 

o 


1.6 


T 




+ cf 






1.2 



(N 



t 0. 





0.4 


+ -H 

o 




T 


1.55 


+ !- 




W 


1.50 






T 
+ ?- 


1.45 


W 


1.40 


+ — 
o 

T 


0.09 


+ ?- 

(N 




pq 


0.08 


1 — . 




T 

+ ?- 

(N 




0.07 







CQ 



0.06 



.a) ' 








' 1 ' 1 




n — 


— □— 


— □— 


— U- 


-n — n— n— n— n-n-n-n-n-n-n-n-n-n-n-n 




■ 








— ■ — p =0 

• p =2 


■ 


- 








-t— p o = 4 
-°— X(5) 


- 


■ 








-*— SU(3) 


■ 








— ■- 

— •- 


■-■■■■■■■■■■■■■■a 


- 


• — 


— •— 


— •— 


| 


* 


=¥= 


=*- 


i i > i 


_ 


b) 








1 i ' i 


- 


■^^ 










■ 


□ — 


^-Lk 


— □— 


— n- 


— n — □— □— □—□-□-□-□-□-□-□-□-□-□-□-□ 


- 






— •— 


-m- 
— •- 


! ! ! ■■■■ 


- 






^— 


— I- 


— I- 


-I- 


i i i i 


- 


c) 








1 i i i 


- 


■^ 


~~-«~_ 








- 


□ — 

• — 


8 


— D — 


— ■- 

— D- 


""■■--■-■-■-■- 

— n — □— □— □— n-n-H = H=n=n=n-n-D.g=a ; g 


■ 


I — 


— I— 


—I- 


I 


- 


i 








i . i 


■ 



10 



15 



20 



Figure 6: (Color online) B(E2) ratios from the 7 band vs. the parameter C as predicted 
by the ES-D model (labelled by the value of the (3q parameter), compared to the X(5) [3], E] 
and SU(3) [29] predictions, as described in subsec. 3.3. 
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Figure 7: (Color online) B(E2) ratios from the /3 band vs. the parameter C as predicted 
by the ES-D model (labelled by the value of the (3q parameter), compared to the X(5) [3], [8] 
and SU(3) [29] predictions, as described in subsec. 3.3. 
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Figure 8: (Color online) B(E2) ratios between the 7 band and (3 band vs. the parameter 
C as predicted by the ES-D model (labelled by the value of the (3q parameter), compared 
to the X(5) [31 [8] predictions, as described in subsec. 3.3. 
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Figure 10: (Color online) (a) Experimental data [45] for the normalized f3 bandhead energies 
Rq/2 = E(0p)/E(2i) vs. the energy ratio i? 4 / 2 = E(4:f)/E(2i), compared to predictions 
of the the ES-D and ES-X(5) solutions, as well as to the predictions of the confined /3-soft 
(CBS) solution [151 EH]- (b) Same as (a), but for the normalized 7 bandheads R 2 /2 = 
E{2~t)/ E{2~i). See subsection 3.5 for further discussion. 
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Figure 11: (Color online) (a) The gamma-stiffness coefficient (3c) 2 is shown as a function 
of the normalized 7 bandhead energy R 2 / 2 = E(2+)/E(2f) predicted by the ES-D solution 
for different values of the Davidson parameter /3q. (b) Same as (a), but for the rare earth 
and actinide nuclei appearing in Table 1 . See subsec. 3.6 for further discussion. 
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Table 1: Comparison of theoretical predictions of the exactly separable Davidson solution 
[ES-D] to experimental data[55] of rare earth and actinides with i? 4 / 2 > 3.0 and known 
0^ and 2+ states. The R4/2 = E(4f)/E(2f) ratios, as well as the (3 and 7 bandheads, 
normalized to the 2^ state and labelled by R 0/2 = E(0^)/E(2f) and R 2/2 = E{2+)/E{2f) 
respectively, are shown. The angular momenta of the highest levels of the ground state, /3 
and 7 bands included in the rms fit are labelled by L g , Lp, and L 7 respectively, while n 
indicates the total number of levels involved in the fit and a is the quality measure of Eq. 
(fTSI) . See subsec. 3.4 for further discussion. 
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Table 1: (continued) 
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Table 2: Comparison of several B(E2) ratios predicted (lower line) by the exactly separable 
Davidson solution [ES-D], for the parameter values shown in Table 1, to experimental 
data[35] (upper line) of several nuclei where the relevant data are known. See subsec. 3.4 
for further discussion. 
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2/3^2i 
2l->0i 

x 10 3 


2fl^4! 
2i->0i 

x 10 3 


2-y >'2\ 

2 7 ^0i 


2 7 ^4i 
2 7 -»0i 


154 Sm 


1.40(5) 

1.47 


1.67(7) 
1.69 


1.83(11) 
1.88 


1.81(11) 
2.06 


5.4(13) 
22.7 


44.5 


25(6) 
142 


1.47 


0.21(5) 
0.08 


156 Gd 


1.41(5) 
1.48 


1.58(6) 
1.73 


1.71(10) 
1.95 


1.68(9) 

2.18 


3.4(3) 
24.6 


18(2) 
52 


22(2) 
179 


1.55(7) 
1.48 


0.16(1) 
0.08 


158 Gd 


1.46(5) 
1.46 


1.68 


1.67(16) 
1.86 


1.72(16) 
2.03 


1.6(2) 
22.1 


0.4(1) 
42.5 


7.0(8) 
133 


1.77(26) 
1.46 


0.079(14) 
0.077 


158 Dy 


1.45(10) 
1.49 


1.86(12) 
1.77 


1.86(38) 
2.02 


1.75(28) 
2.29 


12(3) 
26 


19(4) 
58 


66(16) 
215 


3.22(94) 
1.49 


0.36(15) 
0.08 


160 Dy 


1.46(7) 
1.46 


1.23(7) 
1.67 


1.70(16) 
1.84 


1.69(9) 
2.00 


3.4(4) 
21.3 


40.2 


8.5(10) 
122 


1.89(18) 
1.45 


0.13(1) 
0.08 


162 Dy 


1.45(7) 
1.45 


1.51(10) 
1.65 


1.74(10) 
1.80 


1.76(13) 
1.94 


19.8 


36.2 


104 


1.67(20) 
1.45 


0.14(1) 
0.07 


164 Dy 


1.30(7) 
1.45 


1.56(7) 
1.62 


1.48(9) 
1.75 


1.69(9) 

1.86 


16.9 


29.2 


77 


2.00(30) 
1.44 


0.24(3) 
0.07 


i62 Er 


1.48 


1.74 


1.97 


2.20 


8(7) 
24.9 


52.8 


170(90) 
190 


2.37(25) 
1.48 


0.29(21) 
0.08 


164 Er 


1.18(13) 
1.47 


1.69 


1.57(9) 
1.88 


1.64(11) 
2.07 


22.9 


45.0 


145 


2.19(35) 
1.46 


0.33(5) 
0.08 


166 Er 


1.45(12) 
1.46 


1.62(22) 
1.65 


1.71(25) 
1.81 


1.73(23) 
1.95 


20.1 


36.9 


107 


1.76(18) 
1.45 


0.12(1) 
0.07 


168 Er 


1.54(7) 
1.46 


2.13(16) 
1.68 


1.69(11) 
1.85 


1.46(11) 
2.03 


21.9 


42.1 


131 


1.77(10) 
1.46 


0.129(9) 
0.076 


170 Er 


1.48 


1.73 


1.78(15) 
1.96 


1.54(11) 
2.20 


1.4(1) 
24.8 


0.2(2) 
52.3 


6.8(12) 
184 


1.48 


0.079(19) 
0.080 


166 Yb 


1.43(9) 
1.49 


1.53(10) 
1.77 


1.70(18) 
2.02 


1.61(80) 
2.29 


25.8 


57.4 


215 


1.49 


0.081 


168 Yb 


1.47 


1.70 


1.89 


2.09 


8.6(9) 
23.2 


46.2 


151 


2.09(50) 
1.47 


0.39(10) 
0.08 


170 Yb 


1.47 


1.70 


1.79(16) 
1.90 


1.77(14) 
2.11 


5.4(10) 
23.5 


47.1 


156 


1.78(50) 
1.47 


0.18(5) 
0.08 


172 Yb 


1.42(10) 
1.47 


1.51(14) 
1.69 


1.89(19) 
1.88 


1.77(11) 
2.08 


1.1(1) 
22.9 


3.7(6) 
45.2 


12(1) 
146 


1.48 


0.097(11) 
0.079 


174 Yb 


1.39(7) 
1.45 


1.84(26) 
1.64 


1.93(12) 
1.77 


1.67(12) 
1.89 


18.3 


32.4 


89 


1.45 


0.075 


176 Yb 


1.49(15) 
1.45 


1.63(14) 
1.62 


1.65(28) 
1.75 


1.76(18) 
1.86 


17.0 


29.3 


77 


1.58(11) 
1.44 


0.073 
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Table 2: (continued) 



nucleus 


4i-*2i 
2i-v0i 


6i->4i 
2i-+0i 


8i-»6i 
2i-»0i 


10i-»8i 
2i-»0i 


2i-»0i 

x 10 3 


2i->0i 

x 10 3 


2,3^4! 

2i->0i 

x 10 3 


2-1- — >z\ 
2 7 ^0i 


2 7 ^4 X 
2 7 ^0i 


174 Hf 


1.48 


1.74 


1.96 


2.20 


14(4) 
24.8 


52.5 


9(3) 
185 


1.54(76) 
1.49 


0.0801 


176 Hf 










5.4(11) 




31(6) 








1.47 


1.71 


1.91 


2.11 


23.5 


47.4 


157 


1.48 


0.0791 


178 Hf 


1.47 


1.38(9) 
1.70 


1.49(6) 
1.90 


1.62(7) 
2.10 


0.4(2) 
23.3 


46.7 


2.4(9) 
153 


1.13(17) 
1.47 


0.066(10) 
0.078 


180 Hf 


1.48(20) 
1.48 


1.41(15) 
1.73 


1.61(26) 
1.95 


1.55(10) 
2.17 


24.5 


51.1 


177 


1.34(28) 
1.48 


0.0795 


182 W 


1.43(8) 
1.48 


1.46(16) 
1.72 


1.53(14) 
1.93 


1.48(14) 
2.15 


6.6(6) 
24.1 


4.6(6) 
49.6 


13(1) 
169 


1.98(7) 
1.48 


0.010(1) 
0.0787 


!8 4 W 


1.35(12) 
1.50 


1.54(9) 
1.79 


2.00(18) 
2.07 


2.45(51) 
2.37 


1.8(3) 
26.5 


61.4 


24(3) 
240 


1.91(13) 
1.50 


0.109(9) 
0.0814 


186 W 


1.30(9) 
1.52 


1.69(12) 
1.85 


1.60(12) 
2.18 


1.36(36) 
2.53 


27.3 


69.0 


294 


2.18(15) 
1.51 


0.0829 


186 0s 


1.45(7) 
1.52 


1.99(7) 
1.86 


1.89(11) 
2.19 


2.06(44) 
2.55 


27.3 


70.0 


301 


2.33(12) 
1.51 


0.12(4) 
0.0830 


188 0s 


1.68(11) 
1.53 


1.75(11) 
1.86 


2.04(15) 
2.20 


2.38(32) 
2.56 


27.4 


70.7 


307 


3.20(6) 
1.50 


6.8(13) 
0.0810 


230 Th 


1.36(8) 
1.48 


1.72 


1.94 


2.16 


5.7(26) 
24.3 


50.2 


20(11) 
172 


1.8(8) 
1.48 


0.12(8) 
0.0797 


232 Th 


1.44(15) 
1.46 


1.65(14) 
1.68 


1.73(12) 
1.85 


1.82(15) 
2.03 


14(6) 
21.9 


2.6(13) 
42.0 


17(8) 
130 


2.48(42) 
1.46 


0.045(20) 
0.0770 


234 u 


1.45 


1.64 


1.78 


1.91 


19.0 


34.0 


95 


1.69(40) 
1.45 


0.097(24) 
0.0751 


236 jj 


1.42(11) 
1.45 


1.55(11) 
1.63 


1.59(17) 
1.76 


1.46(17) 
1.88 


17.8 


31.3 


85 


1.45 


0.0743 


238 jj 


1.45 


1.63 


1.45(23) 
1.75 


1.71(22) 
1.87 


1.4(6) 
17.3 


3.6(14) 
30.1 


12(5) 
95 80 


1.74(17) 
1.45 


0.108(12) 
0.0743 


238p u 


1.45 


1.63 


1.75 


1.86 


14(4) 
17.2 


29.8 


11(4) 
79 


1.45 


0.0741 


250 Qf 


1.44 


1.61 


1.72 


1.81 


14.8 


24.7 


61 


1.61(27) 
1.44 


0.092(16) 
0.0730 



34 



Table 3: Values of (3q (determined through minimization of ap n [Eq. (j!5j ) and C [deter- 
mined through Eq. (|14|) ] corresponding to minimum rms deviation G$ a [Eq. (fT5"j) ] between 
the /9i and 71 bands of the ES-D solution for fixed value of R0/2, when the even levels of 
both bands up to L max are taken into account. See subsection 3.7 for further discussion. 



Ro/2 


J-'max 


Po 


C 


a Pn 


5. 


10 


0.00 


1.0 


5.026 


10. 


10 


0.00 


7.3 


3.072 


15. 


10 


1.94 


12.9 


1.141 


20. 


10 


2.64 


16.1 


0.993 


25. 


10 


3.12 


19.4 


0.907 


30. 


10 


3.53 


22.2 


0.855 


5. 


20 


0.00 


1.0 


4.890 


10. 


20 


0.00 


7.3 


4.528 


15. 


20 


1.73 


14.7 


2.611 


20. 


20 


2.57 


17.7 


2.531 


25. 


20 


3.09 


20.6 


2.371 


30. 


20 


3.50 


23.9 


2.197 



Table 4: Theoretical predictions of the ES-D solution for /3o = 1.95 and C = 12.9 compared 
to experimental data for 232 Th [45] . See subsection 3.7 for further discussion. 



L 


gsb 


gsb 


Pi 


Pi 


L 


7i 


7i 




exp 


th 


exp 


th 




exp 


th 





0.000 


0.000 


14.794 


15.021 


2 


15.907 


17.210 


2 


1.000 


1.000 


15.680 


16.021 


3 


16.804 


17.978 


4 


3.284 


3.268 


17.683 


18.288 


4 


18.030 


18.989 


6 


6.749 


6.665 


20.724 


21.686 


5 


19.454 


20.234 


8 


11.280 


11.021 


24.754 


26.042 


6 


21.266 


21.701 


10 


16.751 


16.161 


29.762 


31.182 


7 


23.213 


23.379 


12 


23.033 


21.931 


35.549 


36.952 


8 


25.496 


25.254 


14 


30.035 


28.200 


42.138 


43.221 


9 


27.750 


27.313 


16 






49.438 


49.887 


10 


30.624 


29.542 


18 






57.356 


56.868 


11 


33.219 


31.929 


20 






65.811 


64.102 


12 


36.484 


34.461 
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